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Abstract 

We calculate the bulk viscosity of color-superconducting quark matter in the 
color-flavor-locked (CFL) phase. We assume that the lightest bosons are the su- 
perfluid mode H and the kaons K° and K + , and that there is no kaon condensate. 
We calculate the rate of strangeness-equilibrating processes that convert kaons into 
superfluid modes, and the resultant bulk viscosity. We find that for oscillations 
with a timescale of milliseconds, at temperatures T< 1 MeV, the CFL bulk vis- 
cosity is much less than that of unpaired quark matter, but at higher temperatures 
the bulk viscosity of CFL matter can become larger. 



1 Introduction 



In this paper we calculate the bulk viscosity of matter at high baryon- number density 
(well above nuclear density) and low temperature (of order 10 MeV). On theoretical 
grounds it is expected that, at sufficiently high baryon-number density, three-flavor 
matter will be accurately described as a degenerate Fermi liquid of weakly-interacting 
quarks, with Cooper pairing at the Fermi surface (color superconductivity [1]) in the 
color-flavor-locked (CFL) channel [2]. That is the phase that we will study — the phase 
diagram at lower densities remains uncertain. The motivation for this calculation is 
that in nature the highest baryon-number densities are attained in the cores of compact 
stars, and it is speculated that quark matter, perhaps in the CFL phase, may occur 
there. This means that our best chance of learning about the high-density region of the 
phase diagram of matter is to make some connection between the physical properties of 
the various postulated phases of dense matter and the observable behavior of compact 
stars. Calculating the bulk viscosity of CFL quark matter is part of that enterprise. 

Current observations of compact stars are able to give us measurements of quantities 
such as the mass, approximate size, temperature, spin and spin-down rate of these 
objects. These estimates are steadily improving, and other quantities, such as X-ray 
emission spectra, are becoming available. However, it is a challenge to connect these 
distantly-observable features to properties of the inner core of the star, where quark 
matter is most likely to occur (for reviews, see Ref. [3, 4]). 

One possible connection is via oscillations of the compact star, which on the one hand 
are affected by the transport properties of the interior, and on the other hand may have 
observable effects on the behavior of the star. The bulk viscosity is one of the relevant 
transport properties, and is expected to play an important role in suppressing both 
vibrational and rotational oscillations. One particularly interesting application involves 
r-modes [5, 6, 7, 8]. If the viscosity of the star is too low, unstable r-mode bulk flows 
will take place, which quickly spin the star down, removing its angular momentum as 
gravitational radiation. The fact that we see quickly-spinning compact stars (millisecond 
pulsars) puts limits on the internal viscosity. If we can calculate the viscosity of the 
various phases of quark matter then the observations can be used to rule out some of 
those phases. There are various complications to this simple picture. Viscosity is very 
temperature-dependent, so to obtain useful limits we need good measurements of the 
temperatures of these stars. There are also some uncertainties about additional sources 
of damping that could help to quash r-modes. But the essential point is that viscosity 
calculations of the various phases are of great potential phenomenological importance, 
and in this paper we report on the results of such a calculation. 

The phase that we choose to study is the CFL phase of quark matter. (The bulk 
viscosity for unpaired, non-interacting quark matter has been calculated previously [9, 
10].) It is known that the mass of the strange quark induces a stress on the CFL phase 
that may lead to neutral kaon condensation [11, 12], producing a "CFL-.fr " phase. 
It is not known whether such condensation occurs at phenomenologically interesting 
densities, because of large uncertainties about instanton effects [13], and in this paper 
we will assume that kaon condensation has not occurred: our results are only applicable 
to the CFL phase, where there is a thermal population of K° and other mesons, but no 
condensation. 
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Bulk viscosity arises from a lag in the response of the system to an externally-imposed 
compression-rarefaction cycle. If there are some degrees of freedom that equilibrate on 
the same timescale as the period of the cycle, then the response will be out of phase 
with the applied compression, and work will be done. For astrophysical applications, 
such as r-modes of compact stars, we are interested in periods of order 1 ms, which 
is very long compared to typical timescales for particle interactions. In quark matter 
there is an obvious example of a suitably slowly-equilibrating quantity: flavor. Flavor 
is conserved by the strong and electromagnetic interactions, and only equilibrates via 
weak interactions. 

In unpaired quark matter, the lightest degrees of freedom are the quark excitations 
around the Fermi surface, and their flavor-changing weak interactions produce a bulk 
viscosity [9]. However, in the CFL phase the quark excitations are gapped and their 
contribution to thermodynamic and transport properties at temperatures below the gap 
is irrelevant. Ignoring the quarks, then, the lightest degrees of freedom in CFL quark 
matter are the massless superfluid "if" modes, the electrons and neutrinos, the (rotated) 
photon, and the kaons. Of these, only the kaons carry flavor, so this paper will focus on 
the their contribution to flavor equilibration. In order to calculate the bulk viscosity of 
CFL matter at long timescales such as 1 ms we must therefore calculate the production 
and decay rates of thermal kaons. We expect the dominant modes to be ones that 
involve the H, like K° <-> H H, and ones that involve the leptons, like <-> e v. 

This paper is laid out as follows. Section 2 describes how the bulk viscosity is 
related to the production and decay rates of the kaons. Section 3 describes the basic 
thermodynamics of the system including the kaons and superfluid modes. Section 4 
describes the overall calculation of the rates of the involved processes. Section 5 presents 
the results and conclusions. 

2 Relating bulk viscosity to microscopic processes 

The bulk viscosity is given by [9] 



where the system is being driven through a small-amplitude compression-rarefaction 
cycle with volume amplitude SV (see (2.3) below) and the driving angular frequency is 
to. The average power dissipated per unit volume is 



where r = 2tt/uj. We can parameterize the volume oscillation by an amplitude 5V <C V 
(chosen to be real by convention), and the resultant pressure oscillation p(t) by a complex 
amplitude bp which determines its strength and phase: 



c 



2V 2 dE 
u; 2 (SV) 2 ~dT 



(2.1) 




(2.2) 



V(t) 
p(t) 



V + Re(5Ve iuJt ) 
p + Re(Spe iujt ) 



(2.3) 
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Substituting these into (2.2) and (2.1), we find that 



_ K (2.4) 
Im(<5p) V 

^ ~ 6\^u 

We therefore expect that lm(8p) will turn out to be negative. To determine this quantity, 
we note that the pressure is a function of the temperature and the chemical potentials. 
We assume that heat arising from dissipation is conducted away quickly, so the whole 
calculation is performed at constant T, and in order to find bp we only need to know how 
the chemical potentials vary in response to the driving oscillation. We expect the bulk 
viscosity to be most strongly influenced by the lightest excitations that carry flavor, and 
for the sake of definiteness we will take those to be kaons. Our analysis could easily be 
modified to treat the case where the lightest bosons were pions. At this point we do not 
have to specify whether our kaons are K° or K + . The relevant chemical potentials are 
fid — fi s for the K° and fi u — fi s for the K + . For the following generic analysis we will 
just write the equilibrating chemical potential as "//#■" . 

In thermal equilibrium, the distribution of kaons is determined by their dispersion 
relations (3.13) and the Bose-Einstein distribution. When the kaon population goes 
slightly out of equilibrium in response to the applied perturbation, strong interaction 
processes are still proceeding quickly: only weak interactions are failing to keep up. 
This means that we can always characterize our kaon population by the flavor chemical 
potential fix, so the kaon distribution has the form tik(p) oc p 2 / (exp — px)/T) — 
1) (see Eq. (3.16)), and nonzero [Ik indicates deviation from equilibrium. 

We express the variations in the chemical potentials for quark number and strangeness 
in terms of complex amplitudes <5/i, and 5[1k, 



iwt\ (2-5) 



Note that the term —m 2 /(2fi) which is often described as an "effective chemical poten- 
tial" is already included in the kaon dispersion relation (3.13), so in equilibrium, \ik is 
zero. The pressure amplitude is then 



dp 

OfJ, 



6,+ 9 " 



dfiK 



5[i K = n q 5[i + n K 5[i K , (2.6) 



(From now on all partial derivatives with respect to [i will be assumed to be at constant 
fix, and vice versa.) In principle one might worry that what we have called "n^" is 
really n,i — n s (or n u — n s ), but at temperatures below the gap, and in the absence of kaon 
condensation, thermal kaons make the dominant contribution to — n s and n u — n s . 
From (2.6) and (2.4) we find 

1 V ( _ \ 

C = T77 n q lm(5fi) + n K lm(5fi K ) . (2.7) 

u o V V / 
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To obtain the imaginary parts of the chemical potential amplitudes, we write down the 
rate of change of the corresponding conserved quantities, 



dn q 
~dt 
driK 



dn q dfi dn q dfiK 



dfi dt 
driK dfi 



+ 



d^K dt 
driK dfiK 



dfi dt dfiK dt 



n q dV 
V~dt 
tik dV 



(2.1 



r 



total 



All the partial derivatives are evaluated at equilibrium, fi = p, and fiK = 0. The right 
hand term on the first line expresses the fact that quark number is conserved, so when 
a volume is compressed, the quark density rises. On the second line, which gives the 
rate of change of kaon number, there is such a term from the compression of the existing 
kaon population, but there is also a net kaon annihilation rate of r tota i kaons per unit 
volume per unit time, which reflects the fact that weak interactions will push the kaon 
density towards its equilibrium value. The annihilation rate will be calculated from the 
microscopic physics in section 4. For small deviations from equilibrium we expect r tota j 
to be linear in fiK, so it is convenient to write the rate in terms of an average kaon width 
7#, which is defined in terms of the total rate by writing 



total 



dn K s 
Ik — — Ofi K e 



iuit 



dfi 



A 



(2.9) 



where the derivatives are evaluated at fiK = 0. 

We now substitute the assumed oscillations (2.3) and (2.5) in to (2.8), and solve to 
obtain the amplitudes Sfi and 5 fiK in terms of the amplitude 5V and angular frequency 
uj of the driving oscillation. Inserting their imaginary parts in (2.7) we obtain the bulk 
viscosity 



C = C- 



7cff 



lis 



(2.10) 



where 
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driK (dn 



n K 



n,, 



dfi V dfi 



(2.11) 

The approximate forms on the right hand side are valid for T a. They follow from 
the fact that all the derivatives of the kaon free energy go to zero as T — > 0, so % and 
its derivatives are suppressed relative to n q and dn q /dfi, which are of order fi 3 and fi 2 
respectively. 

To evaluate C and 7 e fj we need the particle densities and their derivatives, which 
follow from the full free energy of the system, Q = f2cFL-quarks(/-0 + Qk((J>, Hk), where 
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^CFL-quarks is the CFL quark free energy at zero temperature [14] and VLk is the kaon 
free energy (3.16). Then n q and dn q /dfi come dominantly from ^cFL-quarks, and all the 
other quantities in (2.11) come from Q^- We can then see that the kaon free energy 
depends on Ek{p) — fiK, and we choose mx to be independent of /x, so from the kaon 
dispersion relation (3.13) we see that the kaon free energy is a function of /i^ + m^/(2/i). 
This means that 

driK_ = _m^_dnK_ 
dfi 2/i 2 dfi K 

so terms in (2.11) involving dnx/d/j, are suppressed relative to those involving dnx/dfix- 
From (2.10) we can already see how the bulk viscosity depends on the angular fre- 
quency uj of the oscillation and the equilibration rate jk- At fixed 7^, the bulk viscosity 
decreases as the oscillation frequency rises; it is roughly constant for uj < 7^, and then 
drops off quickly as 1/uj 2 for to 3> At fixed uj, the bulk viscosity is dominated 
by processes with rate 7^ ~ uj, and their contribution is proportional to 1/7^. If we 
imagine varying the rate but keeping other quantities fixed (e.g. by varying the coupling 
constant of the equilibrating interaction), then for 7^ <C uj or 7^ ^> uj the bulk viscosity 
tends to zero. Thus very fast processes, such as strong interactions, are not an impor- 
tant source of bulk viscosity. The limit of zero equilibration rate and zero frequency is 
singular, and depends on the order of limits. 

In this paper we will also be concerned with the temperature dependence of the bulk 
viscosity. This cannot be straightforwardly read off from (2.10) because the rates and 
particle densities depend on the temperature in complicated ways; however we expect 
that as we go to higher temperatures (T ^> m^, ^k) the bulk viscosity will grow because 
the kaon density is rapidly increasing. In the limit of low temperature we expect the 
viscosity to be suppressed by exp((— nix + /i^)/T) as the thermal kaon population 
disappears. 



3 Dynamics of the light modes 

In this section we lay out the properties of the lightest modes of the system, since they 
will dominate the transport properties. There is an exactly massless scalar Goldstone 
boson associated with spontaneous breaking of baryon number, and some light pseu- 
doscalars associated with the spontaneous breaking of the chiral symmetry. We will 
ignore the 7/ mode associated with the breaking of U(1)a, since U{1)a is explicitly 
broken in QCD at moderate densities. 



3.1 The superfluid "if" mode 



The CFL quark condensate breaks the exact U(1)b baryon number symmetry of the 
QCD Lagrangian, creating a superfluid with an exactly massless Goldstone boson if. 
The Lagrangian for the superfluid mode of the CFL phase is [15] 



L 



N C N 



vli 



12tt 2 



(9 o 0-/i) 2 -(ft0) 2 



(3.1) 



This Lagrangian is correct to leading (zeroth) order in a s and to leading order in the 
derivatives of the field. This can be rescaled to give a conventionally normalized 
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kinetic term, and the total time-derivative term can be dropped [16], giving 



l cS = \{d,<py - \{d4) 2 - ^M^) 2 + ( 3 - 2 ) 

Ignoring the interaction terms for the moment, the dispersion relation for the H particle 
is 

E H (p) = v H p (3.3) 

where v\ — 1/3 is the ratio of the spatial and temporal derivatives above. In thermal 
equilibrium at temperature T, the H bosons have free energy 

T f°° TV 2 



n H = 7T^ dkk ln (! - exp(-v H k/T)) = -t^T 4 (3.4) 
2tt 2 ./n 90vi 



H 



and number density 



d 3 k 1 C(3) , , N 

(In) 6 exp{E H / 1 ) — 1 ■n l v i H 



The cubic and higher-order terms allow a single H to decay into multiple H particles. 
For energies far below [i the dominant process is H — > HH [16, 17], whose rate can be 
calculated by taking the imaginary part of the 1-loop H self energy. Higher order 
corrections to the self energy are ignored, both in Ref. [16] and in this paper. These 
corrections could be calculated by taking into account quantum mechanical interference, 
the Landau- Pomeranchuk-Migdal (LPM) effect. This should only introduce a difference 
of 0(1) in the coefficient of the self-energy but would not change the parametric result. 
The self energy is given by Eq. (3.7) in Ref. [16], 

v < \ 47f2 V" f d " k w( i\( s ^ 2 ^- + f(siE 1 ) + f(s 2 E 2 ) \ 

Zh[Po,P) = -^r^ 1^ / 7T^3 F (Po,P,fc) ttt^ : 4 b , 3.6 

81/x 4 *-^,J \2^Y \4:E 1 E 2 iuj - sxE-i - s 2 E 2 J 



Sl,S2 



(3.7) 



where 

F(po, p, k) = pi — v 2 p 2 — 2vk(p — vk) 

f(E) = \/{e E / T - 1), E l = vk, E 2 = v\p-k\ 

Ref. [16] showed that the real part of this self-energy is parametrically smaller than the 
imaginary part, so we will only concern ourselves with the imaginary part, which we will 
call IL^. There is no contribution to the imaginary part when si = s 2 = —1 as there is 
no pole in the integral. One can also show that the two terms where the signs of Si and 
s 2 are opposite are identical. We can then rewrite this in a slightly simpler and more 
suggestive form that will be used in Section 4. 

tt , \ 2tt 3 p 2 1 f d 3 k 

Uh{P0 ^ = 8W1 + /(P0 ) J 2M^ F ^ P ^ )G ^ P ^ k) ' 

G(p ,p, k) = (1 + f(Eh)) ( 1+ ^ 2) S( Po - El -E 2 )+ 2 ^ 5(p -E 1 + E 2 ] 

(3.1 
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The H propagator can then be written as follows 



D H {p^p) = — r — — — (3.9) 

Po ~ v hP +iTLh{Po,P) 

and we will use this expression in section 4. It will be useful in the calculation of 
the decay rates to have the H self-energy at momenta and energies close to mass shell 
(po = vp). The self-energy is discontinuous at this point so there are two values IT^ and 
11^ depending on whether po tends to vp from above or below, 

n+( P ) = ^n H (vp + e, P ) = - 81 ^ 1 + ) M j o P dki(p,k), 

(3.10) 

2np 1 f°° 
Hi(p) = hrnU H (vp-e,p) = —^—^-^ dkl(p,k). 

where I(p,k) = k 2 (p — k) 2 {l + f{vk))f{yk — vp). As T — > 0, H%(p) oc p 6 //i 4 , and 
Usip) - 0. 

3.2 Pions and Kaons 

The CFL quark condensate breaks the approximate SU(3) chiral symmetry of the QCD 
Lagrangian, creating eight light pseudoscalar pseudo-Goldstone mesons. This octet is 
just a high-density version of the pion/kaon octet. It is described by an effective theory 
[18, 19] 

Los = ^TrfVoXVo^ - + ■ • • (3.11) 

where S = exp(iP a X a /f 7T ), and the normalization of the GellMann matrices is tr(A a A 6 ) = 
25 ab , which yields a conventionally normalized kinetic term for the Goldstone boson fields 
P a . At asymptotic densities, weak-coupling calculations give [19] 

When weak interactions have equilibrated, the pseudoscalars P = tt^, K , K°, K° have 
dispersion relations [19, 11] 



E P = -nf + \jvlp 2 + m 2 P , (3.13) 

where 



^ = 


± m 2 -m 2 u 


2/i 


= 


± m 2 s -m 2 u 


2/i 




m 2 8 - m 2 


2/i 



(3.14) 



Because m s 3> m u , m^, the K and are expected to have the smallest energy gap, 
and so we focus on their contribution to the bulk viscosity. We are interested in study- 
ing small departures from equilibrium, where each meson has an additional chemical 
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potential 5/zp. It will turn out that the K° makes the dominant contribution, so only 
Sfi K o is relevant (Sec. 4.1). 

The expression for the bulk viscosity (2.10) contains terms of the form dnx/dfi, 
which take into account the fact that the meson distributions depend on the meson 
dispersion relations, which via (3.14) depend on the quark chemical potential. In this 
paper we treat the meson masses m^o etc as constants, but in perturbative calculations 
they also depend on \i and the CFL pairing gap A (see section 3.3). 

In thermal equilibrium at temperature T and with chemical potential /ip, the free 
energy and number density of a meson P is 



T 

tt P -- 



2ir 2 



POO 

/ dkk 2 ln(l-exp(-(£p-5/ip)/T)) (3.15) 
Jo 



oo 



n P = — rs = — k / dk k 2 — — — — (3.16) 

When weak interactions have equilibrated 5fip = 0, but when weak interactions are out 
of equilibrium the mesons may have nonzero chemical potentials. 



3.3 Pseudo-Goldstone-boson masses 

Although we treat the masses as constants, they are predicted to have density depen- 
dence in high density QCD [13] 

ml ± = -fi(2A + ABm s ) {m u + m d ) , 

J n 

m 2 K ± = -^{2A + 4Bm d ){m u + m s ) , (3.17) 

J n 

where A is positive and related to instantons [20, 21]. In the limit of asymptotically 
large density the coefficient A can be computed reliably, but at moderate density its 
value is quite uncertain [13]. Asymptotic-density QCD calculations [19, 11] also yield 

3A 2 

although it is not clear how well these expressions can be trusted at densities of phe- 
nomeno logical interest. We will assume that A and B are such that there is no meson 
condensation at zero temperature, which means that all the meson masses are greater 
than their effective chemical potentials. 



3.4 Weak interactions between light bosons 

We have argued above that the bulk viscosity will arise from flavor violation, which 
will be dominated by conversion between the lightest pseudo-Goldstone modes (neutral 
kaons, typically), which carry flavor, and the superfluid H modes, which are flavorless. 
The dominant effect of the weak interaction will be to introduce mixing between the 



8 



K° and the H, in the form of a K° <-> vertex (3.26) in the effective theory. We now 
calculate the strength of that coupling. 

The Lagrangian density for the H modes can be written in a nonlinear form analogous 
to (3.11) for the pseudoscalars, so the leading terms in the CFL effective theory become 
[18], 

£ eff = ^Tr [VqXVo^ - ^Ec^] + 12f 2 H [V ZV Z* - v^Zd^*] (3.19) 

where Z = exp(iH/2\/6fH) is the field related to the breaking of U(1)b and fn is the 
corresponding decay constant. At large density the coefficients of the CFL effective 
Lagrangian can be determined in perturbation theory. At leading order v 2 . = v 2 H = 1/3, 
f n is given by (3.12), and from Ref. [19] we obtain 



4 \2tx\ 

Under an element (L, R, expect) ) of the the chiral flavor and baryon number sym- 
metry group SU(3)l x SU(3)r x U(l)v, the left-handed quarks, right-handed quarks, 
and bosons transform as follows: 

q L — > exp(ia)Lq L , 

q R -> exp(ia)Rq R , . , 

Z — > exp(ia) Z . 

The weak Hamiltonian breaks the approximate flavor symmetry of QCD, and only acts 
on left-handed fields. The elementary process that is relevant to kaon decay is the 
conversion between strange quarks and down quarks via exchange of a W^, which can 
be treated at energy scales well below 100 GeV as a four-fermion interaction (Ref. [22], 
sect. II-3 and II-4), 

£ weak = GpV ^ us (l7 M n) L (n7 M c?) L + h.c. (3.22) 

where V u dV us ~ 0.215(3). In order to determine how this interaction is represented in the 
low energy effective theory of the CFL phase, we introduce the spurion field A^ s which 
transforms as A^ s — > LA^LL In the QCD vacuum we set A^ s = A6 using the usual 
notation for the GellMann matrices (Ref. [22], sect. II-2), i.e. (A^)^ = 62083/3 + d-sa^p- 
Thus each time this spurion field occurs in an interaction, it mediates a conversion of 
downness into strangeness, or vice versa. The lowest-order terms in the effective theory 
that involve such a conversion are obtained by writing down the lowest-order terms that 
contain A ds and are invariant under spatial rotations and SU(3)l x SU(3)r x U(1)v- 

£wcak = flf H G ds Ti [A ds (sa E t Zd Z* - vimtf ZdiZ*^\ (3.23) 

where Gd s and Vd s are new couplings in the effective action. 

Dimensional analysis suggests that Gds ~ Gf and Vd. s ~ 1. If the density is large 
we can be more precise and determine the coupling constants using a simple matching 
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d,s u 




s,d u 



Figure 1: Leading contribution to the X®B U polarization function in the microscopic 
theory with gauged chiral and baryon number symmetries. The shaded bar corresponds 
to the vertex of (3.22), which is the low energy limit of a VF-boson-mediated interaction. 



argument. For this purpose we gauge the SU(3)l and U(1)b symmetries. We will 
denote the corresponding gauge fields and B^. The flavor-violating term (3.23) in 
the effective action leads to mixing between the X® and gauge bosons. By matching 
to a calculation of the mixing term in the microscopic theory (3.22), we now proceed to 
determine G ds in terms of Gp. 

In the effective theory, the E field has one left-handed quark index so c^E — > (d^ + 
X^ Aa)S. For the superfluid mode, d^Z — > (<9 M + B^Z. Substituting into (3.23) and 
evaluating in the CFL vacuum (E = 1) we find the mixing term is 

C mix = 2G ds f 2 f 2 H (X«B - viXfB,) . (3.24) 

In the microscopic theory, the corresponding calculation is the computation of the X®B U 
polarization function. At weak coupling the dominant contribution comes from the two- 
loop diagram shown in Fig. 1. The evaluation of the Feynman diagram is described in 
Appendix A. The result (A. 5) is 

G ds = V2V ud V us G F v 2 ds = v 2 = 1/3. (3.25) 

It is now straightforward to read off the K° — > H amplitude. Linearizing (3.23), we find 

C = G ds fJ H {d K°d H - vjA^diH) (3.26) 

with G ds and v ds given in (3.25). This leads to a vertex factor for the K-H interaction 
given by 

A = G ds f 7T f H (pl-v ds p 2 ); (3.27) 

This is the value of the K°-H vertex in Feynman diagrams such as Fig. 2. Combining 
this vertex factor and the Lagrangian for the H, we can calculate the matrix element 
for conversion between a kaon with 4-momentum p and two Hs with 4-momenta k and 

M 2 K o HH (P, k, q) = G2Jl { fl~ f f dsP2? (po(k ■ q) + k {p ■ q) + q {p ■ ktf \D H (p ,p)\ 2 , 

(3.28) 

where D H is the if -propagator (3.9). 
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H 
H 



Figure 2: The K° — > H H diagram, including the K° — > H vertex (square), a full 
H propagator (thick line) and the H — ► H H vertex (round). All external lines are 
amputated. 



4 Rates of strangeness re-equilibration processes 
4.1 Neutral Kaon Rates 

In principle, the correct way to calculate the K° annihilation rate is as follows. As noted 
above, the weak interaction introduces a small mixing between the K° and the H. We 
rediagonalize the kinetic terms in the effective action, in terms of new fields Ek (which 
is the kaon with a tiny admixture of H) and Eh (which is the H with a tiny admixture 
of K°). We have already seen that the H has a width, arising from the possibility of 
H — > H H 1 . When we diagonalize, this will induce a width for the Ek- Because the 
Ek is almost the same state as the kaon, the Ek width is a very good estimate of the 
K° width. In terms of the original basis, this width arises from the vertex shown in 
Fig. 2. In the interests of brevity, we do not rediagonalize, but simply calculate the 
contribution of the vertex shown in Fig. 2 to the kaon annihilation rate. This comes via 
the processes K° «-> H H and if K° <-» if. The net kaon annihilation rate is 

r — r _ r — n _ *>-<Wo/ T, ir ~ (a-\\ 

1 total 1 forward 1 backward l, 1 C JL forward ~ 7 1 forward > v^ - - 1 -/ 

where we have used the properties of the Bose-Einstein distributions. We keep only first 
order in 5{j,k°, and obtain the average kaon width 7^ (2.9), remembering that 5fiK in 
section 2 was the amplitude of a complex oscillation, so S/ik = SfiKexp(iut), 



( dn K \ 1 T {orward (5fi K = 0) , 4 2 j 

\dji K ) T 



We can therefore obtain the average kaon width simply from the forward rates K° — ► 
if and if if — > ii. The contribution from if — > H H is 

rxo^H = \ff! \ M \ 2 ( 27r ) 3 5 (P - 91 - 92) (2tt) 5(po - vgi - vq 2 ) F BE (p , q u q 2 ) 

Jp J qi J (J2 

(4.3) 

with |M| 2 given in (3.28) and F BE (p , q u q 2 ) = f(p - Sfi K o)(l + /(ugi))(l + f(vq 2 )). 
The rate for i< ii — ► if can be obtained by multiplying by 2 for the symmetry factor 
difference, switching q 2 — > —q 2 in both delta functions and turning (1 + f VQ2 ) f vq , 2 to 



1 There are also decays involving three or more H particles, but we expect these to be suppressed, 
since the H is derivatively coupled, and the greater the number of H particles involved, the smaller the 
momentum carried by each of them. 
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make that H an incoming particle. Adding the two contributions, and performing the 
q2 integral using the momentum-conserving delta-function, we obtain the total forward 
rate 



forward 



GMifH f(E K )(l + f(E K )) 



{El-v±p 2 fXl H {E K ,p) 



(E 2 K -v^ P 2 y + u H (E K , P y 

where Hh(Ek,p) was defined in Eq. (3.8). 



Forward Rate Integrand 
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(4.4) 



Figure 3: Plot of the integrand of (4.4) as a function of momentum for p < p. In 
this plot fi = 400 MeV, m s = 120 MeV, and m K = 27.92 MeV, so p = 22 Ah MeV 
and T a ~ 0.14 MeV. For T ^ p the integral is dominated by low-momentum kaons, 
i.e. p p (T = 0.05 MeV line in plot). But as T gets closer to p, the near-singularity at 
p = p becomes more important, and at T = 0.5 MeV the integral is actually dominated 
by the region where p is very close to p. 



In general this can be evaluated numerically and then combined with (4. 2), (2. 10), (2. 11) 
to obtain the bulk viscosity. However for certain temperatures, the integrand is dom- 
inated by momentum that make the denominator as small as possible, i.e. when 
Er = vp. This corresponds to the H «-> K° resonance, at which the kaon has a 
special momentum p such that an H with that momentum is also on shell, 



= rg* = _ i + . (4.5) 

where 5m = rriK — fJ^o- Near this momentum, the virtual H in Fig. 2 is almost on 
shell, so momenta close to p dominate the integral. As long as the numerator is slowly 



12 



varying, we can approximate the sharp peak as a delta-function, obtaining 



g 2 f 2 f 2 9 p v p/ t 

f ™ d * (1 + </^ V (e ^_ 1)a • (4-6) 

This expression becomes invalid at very low temperature T<cT a , when there are very 
few kaons with momentum p so the main contribution does not come from p pa p, and 
at very high temperature T ^> T& when there are so many thermal kaons with p > p 
that they outweigh the contribution from the resonance. One determines T a and T& 
by setting the non-resonant contribution equal to the resonant value from Eq. (4.6), 
giving T a7 b as a function of 5m. Tb has to be determined numerically and we found that 
Tb ~ 9 MeV for 5m = 0.1 MeV and monotonically increases as 5m increases, so Tb is 
almost always higher than the temperature range that is of physical interest. T a is given 
by the following condition 

Sm2 , (4.7) 



(K 



a 2pL% V \6m{m K T a )W 

where the T a dependence on the right side is logarithmically weak, so we expect that 
T a <5m 2 /(2^ ). 

The integrand of (4.4) is plotted in Fig. 3, for several different values of the temper- 
ature, showing that for lower temperatures, T < T a , the integrand is a smooth function, 
with a broad peak in the low-momentum region, but as the temperature rises, and the 
number of thermal kaon with momentum p rises, the integrand develops a very sharp 
peak at p — p, where the intermediate H is on shell. 



4.2 Charged Kaon Rates 

In principle there is also a contribution to kaon number violation from the charged kaon 
modes, which necessarily involve charged leptons. We will now show that this can be 
neglected compared to the contribution from the neutral kaons. The lightest charged 
kaon is the K + , and the relevant creation/annihilation reactions are 

K+ «-> e + u e (4.8) 
K + + e~ w v e (4.9) 
K + + v e <-> e + . (4.10) 

The matrix element for this process has been calculated in Ref. [23, 24] and is given by 

A = GUsmO^eihh^l ~ 75)^2), (4.11) 

where G is the appropriate coupling constant for the charged kaons in the medium, 
which we expect is of order Gp. Summing over initial spins and averaging over final 
spins, we find 

M 2 = G 2 fl sin 2 6 C m 2 e (h ■ k 2 ) (4.12) 
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The rate of the first reaction is 



T = r dx r dyxy X + V + ^ +/T (\ - z )F(x, y) (4.13) 



o J yi 



x + y 



F(x,y) 



1 

2v 2 xy 
8fi K+ e x+y 



z = ^ Z7 \{x 2 + y 2 ){l-v 2 ) + 2{x + y)iif + /T + 2xy + {iif + 2 -m 2 K+ )/T 2 



T (e x +y - l){e x + l)(ey + 1) 



where F(x, y) is the product of the distribution functions for the kaon and two leptons to 
lowest order in Sfix+- The rates for the second and third reactions, which are identical 
for a massless electron, can be derived from a simple change of x — > — x and y — > —y, 
respectively. They can then be evaluated numerically. Note that these calculations are 
done keeping only the lowest order term in m e and setting \i e — 0. 

We can then compare this rate to the rate for neutral kaon decay and find that 

Tk+ ~ -. (4.14) 



2 



for m>K0 ~ rn>K+, so the contribution from charged kaons is suppressed by a factor of 
(T//x) 2 . This is to be expected, since the phase space for quarks is of order /i 2 T, localized 
near the quark Fermi surface, whereas the phase space for electrons is of order T 3 , since 
in the CFL phase there is no Fermi sea of electrons. 



5 Results 

Our result for contribution of kaons to the bulk viscosity of CFL quark matter is given 
by equations (2.10), (2.11), (4.2), (4.4). The bulk viscosity is most sensitive to the 
temperature, and to the kaon energy gap 

2 2 

Sm = m K - fi e £ = m K s — . (5.1) 

It also depends on the orthogonal combination rrix + /^o, the quark chemical potential 
/i and CFL pairing gap A. As discussed in section 3.3, we have no reliable way to 
calculate mx in the density range of interest for compact stars, so in presenting our 
results we will treat Sm and T as parameters. 

Using the definitions of the densities of kaons and quarks, one can derive the asymp- 
totic versions of C for temperatures far above and far below the kaon energy gap (Table 
1). One can also derive the low temperature version of the rate and from that, combined 
with C, we can get the low temperature version of the bulk viscosity (Table 2). As one 
would expect, most kaon-related quantities, including the bulk viscosity, are suppressed 
by exp(— 5m/T) at low temperatures. This is because the energy gap Sm is the mini- 
mum energy required to create a K°, so the population of thermal kaons is suppressed 
by a Boltzmann factor. 

To illustrate the likely contribution of kaons to the bulk viscosity of quark matter in 
compact stars, we now evaluate the bulk viscosity numerically for a range of Sm and T. 
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quantity 


asymptotic form 




T < 5m 


m K «T«/i 


n K 


{m K Tfl 2 e- Sm l T 


T 3 


dn K 
^0 


m K (m K Tf/ 2 e- 5m / T 


T 2 


dn K 




m K rp2 


n q 






dn q 






c 







Table 1: Asymptotic forms for the densities and the C parameter. Constant numerical 
factors are not shown and it is implicitly assumed that T < 0.57A so that there is a 
CFL condensate, even when T 3> m K . 



quantity 


T < T a (5m) < 5m 


approximate form 

T a (5m) <T <5m 


r forward 


G\ ^m 3 K T 3 Sm 5 e- dm /' r 






7eff 


G\ 5m 5 




G 2 F // p 4 (m K T)" 3 / 2 e - 5m2 /( 2 <o)T 


c 


Gldrrfm^T^e- 5 ™' 7 uj- 


-2 


G 2 F fiY m 2 K T" 1 e-^ T /( 7c 2 ff + u 2 ) 



Table 2: Approximate forms of the bulk viscosity and related quantities, for small T. 
Constant numerical factors are not shown. The rate has two separate ranges within the 
T <C 5m region: T < T a <C 5m and T a < T <C 5m, where T a < 5m 2 /(2^ ) (4.7). 
Note that p is related to 5m by (4.5). The low temperature entry for ( is in general 
proportional to (7 2 ff + u 2 )~ l rather than just u~ 2 , but in this temperature range 7 c fj 
is always much less than astrophysically relevant frequencies (ou > 1 Hz). There is no 
third column for the higher end of the range of temperatures that we study in this paper, 
T ~ mx, /i^o, because although we can still use (4.6) for rf orwar d, there is no simple 
form for and hence for 7 er r or £. 
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Our calculations are performed at fi = 400 MeV. We vary 5m by varying with fj?^ 
fixed at 17.92 MeV, corresponding to m s = 120 MeV. 

Compact stars have internal temperatures in the MeV range immediately after the 
supernova, and then cool to temperatures in the keV range over millennia, so we explore 
the range 0.01 MeV < T < 10 MeV. Since mx and are both expected to be of 
order tens of MeV [13], we expect 5m to be generically of the same order, so we explore 
the range 0.1 MeV < 5m < 10 MeV. 

The bulk viscosity is determined by the kaon equilibration rate 7^ (4.2) and the 
coefficient C (2.11), so we plot these quantities separately before plotting the bulk 
viscosity. 

5.1 Proportionality constant C (Fig. 4) 
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Figure 4: Coefficient C (2.11) as a function of 5m (left panel) and temperature (right 
panel) 

In Fig. 4 we show how C depends on 5m and T. Roughly speaking, C measures how 
sensitive the kaon and quark number are to changes in \ik and \i. At low temperatures 
T <^ 5m, C is suppressed by an exponential factor exp(— 5m/T), so the curves drop 
rapidly in the high 5m region of the left panel, and the low T region of the right panel. 
We also see that the curves for different 5m start to converge at high temperature (right 
panel). This is because at high enough temperature (beyond the range that we study) 
C would become proportional to T 4 , independent of 5m (see table 1). 



5.2 Kaon width 7 cff (Fig. 5) 

In Fig. 5 we show how the neutral kaon effective width 7^ depends on T and 5m. (We 
also show one charged kaon width curve to illustrate that it is subleading (4.14)). 

The 5m- dependence is shown in the left panel. From Table 2, we expect that at 
a fixed temperature T, for sufficiently large 5m, T a (5m) will become greater than T, 
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Figure 5: Plot of average K° decay width 7^ (4.2), (4.4) as a function of 5m (left panel) 
and temperature (right panel). The horizontal dashed line shows where the width is 1 
kHz (u/2ir = 1 ms^ 1 ), the fastest rotation rate of compact stars. The charged kaon 
width is also shown (dotted line) to illustrate that it is a subleading contribution to 
strangeness equilibration (4.14). The transition that occurs at 5m rs T is where the rate 
becomes dominated by the H resonance (Section 4.1). 



and 7x will then rise as 5m 5 . This is seen at the upper end of the T = 0.01 MeV 
curve, and corresponds to the region where low-momentum (p < p) kaons dominate 
the rate. For the rest of the T = 0.01 MeV curve, and for all the other curves in the 
plot, the equilibration is dominated by kaons at the if- resonance, with momentum p. 
The width shows a peak as a function of 5m, which follows from the approximate form 
for 7 c g given in table 2 (second column). Using (4.5) to relate 5m to p, one can see 

7efr ~ <5m 4 exp(— 5 m 2 /(2/j c ^ Q T)), which is peaked at 5m pcak = 2^J ^ T. For our plots 

H™* m 30 MeV, which gives the observed positions of the peaks. 

The T-dependence is shown in the right panel. At the very lowest temperatures 7 e g 
will have a constant value which depends on 5m (Table 2). This is clear in the curves 
for 5m = 5 MeV and 5m = 1 MeV. As with the large 5m region of the left panel, 
this is where the low momentum kaons are dominating the rate. In the intermediate 
temperature region the width rises quickly and then peaks and drops off slowly. This 
comes from competition (4.2) between Ff orwar d/T, which is monotonically increasing with 
T (4.6), and {dnx/ ' d[iK)~ x , which is monotonically decreasing with T. At high enough 
temperature, The expression (4.6) for Ff orwar d/T rises as T, while drix/dfiK rises more 
quickly, so the width drops. 

At high enough T, the curve with 5m = 0.1 MeV starts to bend upwards. This 
feature actually corresponds to T > T& (from Section 4.1) where (4.6) becomes invalid 
and kaons with high momentum (p > p) dominate the rate. For the other curves, Tf, is 
beyond the range that we study. 
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5.3 Bulk viscosity ( (Fig. 6 and 7) 
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Figure 6: Plot of bulk viscosity as a function of 5m (left panel) and temperature (right 
panel) . 

It is useful to compare the behavior of the bulk viscosity in CFL quark matter with 
its behavior in quark matter phases whose re-equilibration is dominated by ungapped 
fermionic modes such as the 2SC or single- flavor phases (see solid black line in Fig. 6). 
In such phases, the bulk viscosity shows a peak as a function of T, because 7 e g varies 
monotonically with T, while C is determined by the phase space at the Fermi surface, 
and hence is insensitive to T [9, 25, 26]. This produces a single peak when 7 e g- is 
equal to the angular frequency u> of the applied compression oscillation (see (2.10)). As 
we will describe below, our results show that in the CFL phase, the situation is more 
complicated, because 7 c g is no longer a monotonic function of T, and also because C 
can vary rapidly as the control parameters 5m and T are varied. 

The dependence of the bulk viscosity (at frequency uj/2-k = 1 kHz) on the kaon 
energy gap 5m is shown in the left panel of Fig. 6. From consideration of the factor 
of 7eff/(7eff + ^ 2 ) i n (2.10) we would have expected two peaks for T = 0.01 MeV and 
0.1 MeV, because from Fig. 5 we see that at these temperatures 7 e ff passes through 
uo=l kHz at two different values of 5m. In fact we get one peak, close to the lower 
value of 5m at which 7 c fj = u>. The higher peak is washed out by rapid variation of C 
with 5m, which occurs when 5m > T (see Fig. 4). Even outside the physically relevant 
range of 5m shown in our plots, we do not find additional peaks in the bulk viscosity. 

The dependence of the bulk viscosity (again at uj/2tt = 1 kHz) on the temperature 
T is shown in the right panel of Fig. 6. It is a monotonically increasing function of 
T for all values of 5m. This is because, as is clear from the right panel of Fig. 4, C 
varies rapidly with temperature for all physically relevant values of 5m. In fact, the 
temperature-dependence of C dominates the bulk viscosity, so the right panel of Fig. 6 
looks similar to the right panel of Fig. 4. 

Finally, Fig. 7 shows a plot of bulk viscosity ( as a function of temperature for 
different oscillation timescales, r = 2-k/uj. We see that for unpaired quark matter, 
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Figure 7: Bulk viscosity as a function of temperature for oscillations of different fre- 
quencies. The curves peaked on the left are for unpaired 3- flavor quark matter [9] and 
the rising curves on the right are our calculation of the kaonic bulk viscosity of CFL 
quark matter, for 5m = 1 MeV. 

Cunp is independent of to at high temperatures, because 7 unp then rises far above lo, 
so, by (see (2.10)), ( = C/7 unp . However, (cfl depends more strongly on r at high 
temperatures, because 7 e fj is not much greater than uj at high temperature (Fig. 5), so 
the a>-dependence in (2.10) is not suppressed. The CFL bulk viscosity becomes larger 
as the frequency drops. 

6 Conclusions 

We have calculated the contribution of the lightest pseudo-Goldstone bosons, the neutral 
kaons, to the bulk viscosity of CFL quark matter. Our results are given by equations 
(2.10), (2.11), (4.2), (4.4), and are displayed for reasonable parameter choices in Figs. 6 
and 7. The bulk viscosity is most sensitive to the temperature, and to the kaon energy 
gap 5m (5.1). We find that, as one would expect, the kaonic bulk viscosity falls rapidly 
when the temperature drops below the kaon energy gap, since the the kaon population 
is then heavily Boltzmann-suppressed. It is clear from the right hand panel of Fig. 6 
that once the temperature falls below the 10 MeV range (which is expected to occur in 
the first minutes after the supernova [27]) the bulk viscosity of CFL quark matter at 
kHz frequencies is suppressed by many orders of magnitude relative to that of unpaired 
quark matter. 

It is noticeable that at low temperatures the suppression is less severe for smaller 
kaon energy gaps. However, 5m is a poorly-known parameter of the effective theory 
of the pseudo-Goldstone bosons. It is the difference of the kaon mass and the kaon 
effective chemical potential, both of which are expected to be roughly of the order of 
10 MeV [13], so it is unnatural to assume that 5m is much smaller than an MeV or 
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so. For astrophysical applications, it is clear that CFL quark matter can be sharply 
distinguished from quark matter by its bulk viscosity (as by many other transport 
properties) after the very earliest times in the life of a compact star. Also, a rapidly 
vanishing bulk viscosity as the temperature drops below 10 MeV could be a potential 
observable associated with core collapse supernovae. 

In general, bulk viscosity arises from re-equilibration in response to compression. 
We have calculated the dominant contribution to the re-equilibration of flavor in quark 
matter, and we believe that this is the dominant contribution to the bulk viscosity as 
a whole in the range of frequencies that are of astrophysical interest, namely zero to 
1000 Hz. Any other contribution would have to come from degrees of freedom that 
equilibrate on a similar timescale, and the only possibility that we can imagine is the 
thermalization of the low-momentum tail of the thermal distribution of H particles. 

Our results highlight several interesting questions for future research. A natural 
next step would be to extend our calculation to the kaon-condensed "CFL-i^°" phase, 
which corresponds to allowing the kaon energy gap to drop to zero. There are also some 
technical issues in our calculation of the flavor changing rate that remain to be addressed. 
The graph shown in Fig. 2 includes the width of the H boson due to the one-loop thermal 
self energy. This corresponds to the resummation of a class of diagrams with multiple 
H boson radiation and absorption [16]. However, since the mean free path associated 
with small angle two-body collisions is of the same order of magnitude as the radiation 
length (the mean free path between if-bremsstrahlung events) this approximation is 
not correct. A more complete approach has to take into account quantum mechanical 
inteference, the Landau-Pomeranchuk-Migdal (LPM) effect, between different diagrams 
that have the same final state [28]. Another relevant improvement would be to include 
higher-order corrections to the H dispersion relation [29] which could have a strong 
effect on the collinear splitting amplitude for H particles. We do not expect that these 
improvements will affect our results significantly, but they would change the numerical 
prefactors in the rate. Finally, it should be noted that we treated the kaon mass as a 
numerical parameter, but it is expected to be density- dependent, and its //-dependence 
will feed into our expressions for the bulk viscosity. We expect that this would only 
weakly affect our results, but we have not performed an explicit check. 
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Appendix A Matching calculation of Gds and 



In this appendix we explain the matching calculation of the coefficients G^ and v ds in 
the effective lagrangian. As discussed in Sect. 3 these two coefficients can be extracted 
from the — X® polarization function at zero energy and momentum transfer. The 
leading order contribution to the polarization function at very large baryon density is 
shown in Fig. 1. This two- loop diagram in can be viewed as the product of two one- loop 
polarization functions. This can be made more manifest by Fierz rearranging the weak 
interaction given in (3.22), 

C = - GFV ^ Vus (sj,d) L (urn) L . (A.l) 
We can now write the B^-X® polarization function at uj = q = as 

r xb = -^«fn^nw^, (A.2) 

where Tlffi L and U.™ are the polarization functions of the currents (sj^d^L and 
(u / ~f fl u) l in the limit of vanishing energy and momentum transfer. The polarization 
functions have the tensor decomposition 

IV = f (5 fl0 5 u0 - v 2 5 fli 5 l/i ) , (A. 3) 

where / and v can be viewed as the decay constant and velocity of a collective mode. 
The relevant correlation functions in the CFL phase were first computed by Son and 
Stephanov [19]. Using their results we can write 

K B = -^^p^VlPn (W*> - A,**) , (A.4) 

where f w and fjj are the pion and H decay constants, see equ. (3. 12), (3. 20), and v is the 
Goldstone boson velocity. Comparing with (3.24) we conclude that 

G ds = -V2G F V ud V us , v 2 ds = v 2 = 1/3. (A.5) 
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